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ABSTRACT 

We suggest a possible explanation for the high frequency quasi-periodic os- 
cillations (QPOs) in black hole low mass X-ray binaries. By solving the pertur- 
bation general relativistic magnetohydrodynamic equations, we find two stable 
modes of the Alfven wave in the the accretion disks with toroidal magnetic fields. 
We suggest that these two modes may lead to the double high frequency QPOs 
if they are produced in the transition region between the inner advection dom- 
inated accretion flow and the outer thin disk. This model naturally accounts 
for the 3 : 2 relation for the upper and lower frequencies of the QPOs, and the 
relation between the black hole mass and QPO frequency. 

Subject headings: magnetohydrodynamics - QPOs - accretion disc - stars: black 
hole 

1. Introduction 

Low-mass X-ray binaries (LMXBs) are binary systems consisting of a neutron star (NS) 
or black hole (BH) accreting from a low-mass (< 1M Q ) companion star. X-ray emission of 
LMXBs often shows fast X-ray variability in the form of high frequency quasi-periodic oscil- 
lations (HFQPOs), which frequently appear in pairs in certain state simultaneously (van der 
Klis 2006). Abramowicz & Kluzniak (2001) pointed out that the frequency ratio of the twin- 
peak HFQPOs in the BH source GRO J1655-40 equals 3/2, and that this commensurability 
of frequencies may be a signature of a non-linear resonance. Later, Abramowicz et al. (2003) 
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found a signature of the same commensurable ratio in the twin-peak HFQPOs observed in 
an NS-LMXB, Sco X-l. Based on this observational evidence, Kluzniak and Abramowicz 
argued in several papers (Kluzniak & Abramowicz 2001; Kluzniak, Abramowicz, & Lee 2004; 
Kluzniak et al. 2004) that the twin-peak HFQPOs in both BH and NS sources are due to 
the same physical mechanism — a non-linear parametric resonance in accretion disk global 
oscillations. 

However, while for the BH sources the commensurable ratio 3/2 was quickly confirmed 
and generally accepted (e.g. Remillard & McClintock 2006), the presence of the same com- 
mensurability in the NS sources is denied by several experts (e.g. Boutelier et al., 2009). 
There is no consensus whether the nature of the twin-peak HFQPOs in the two types of 
LMXBs is the same. We have proposed a mechanism for the twin kilohertz QPOs in NS- 
LMXBs using the magnetohydrodynamic (MHD) Alfven wave oscillations, and the results 
seem to fit the observation well (Li & Zhang 2005; Shi & Li 2009). In this paper we focus 
on an MHD explanation of the HFQPOs in BH-LMXBs. 

Barret et al. (2005) measured the quality factor Q for the HFQPOs measured in the 
NS-LMXB 4U 1608-52, and found that Q ~ 200. They argued that such high coherency 
is impossible to achieve from kinematic effects in orbital motion of hot spots, clumps or 
other similar features located at the accretion disk surface, because these features are too 
quickly sheared out by the differential rotation. Although orbital motion cannot explain the 
the HFQPOs in LMXBs, the frequencies of several fluid oscillatory modes are expressed by 
the three characteristic orbital frequencies: the "Keplerian" frequency, the "radial" epicyclic 
frequency, and the "vertical" epicyclic frequency. In the Kerr metric, these three orbital 
frequencies and the Lense-Thirring "frame-dragging" frequency have been listed (e.g. Perez 
et al. 1997). Several HFQPOs models use their ratios (in various combinations) to explain the 
observed 3/2 commensurability. Cui, Zhang, & Chen (1998) suggested the Lense-Thirring 
nodal precession frequency near the inner stable circular orbit (ISCO) radius as the lower 
HFQPO frequency, such as the 300 Hz QPOs in GRO J1655-40. The relativistic precession 
model of Stella et al. (1999) applies to both BH and NS sources; the pariastron precession 
frequency and the Keplerian frequency were taken as the lower and upper frequencies of the 
twin HFQPOs, respectively, whereas the QPOs at much lower frequencies were interpreted in 
terms of the Lense-Thirring nodal precession ferquency. Wang et al. (2003, 2005) suggested 
that a non-axisymmetric magnetic coupling of a rotating BH with its surrounding accretion 
disk coexists with the Blandford-Znajek process. The two frequencies were supposed as the 
Keplerian frequencies of two hotspots, one near the inner edge of the disk and the other 
somewhere outside, respectively. 

Wagoner et al. (2001) considered the modes of the diskoseismic wave, such as g-modes, 
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p-modes and c-modes as the explanation of the HFQPOs. They estimated the masses and 
angular momenta of some BHs with the measured frequencies of the HFQPOs when the g- 
modes or c-modes were selected. Rezzolla et al. (2003) discussed the inertial-acoustic modes 
in a small-size torus very close to the horizon of the BH while the centrifugal and pressure 
gradients were selected as the only restoring forces. In this model the black hole spin had 
to be very close to the maximal value to produce the 3 : 2 ratio. Tassive & Bertschinger 
(2008) investigated the kinematic density waves in the accretion disks when nothing but the 
gravity was considered as the restoring force, and several discrete radii were adopted. Several 
modes in pairs close the ratio (3 : 2) could be got, but the correct frequencies couldn't be 
reproduced. Similar to the parametric resonance models, it is difficult to explain why other 
modes such as the fundamental frequencies weren't observed except the two modes in pairs. 

In modeling HFQPOs in BH-LMXBs two points need to be mentioned. First, in most 
models the 3 : 2 ratio was often overemphasized and substituted into those models directly. 
In fact, the 3 : 2 ratio of the twin HFQPOs in BH-LMXBs isn't rigorous but approximate. 
Second, the HFQPO frequencies were often considered invariable so that these frequencies 
(168, 113 Hz and 67, 41Hz in GRS 1915+105, Remillard (2004); 450, 300 Hz in GRO J1655- 
40, Strohmayer (2001a); 276, 184 Hz in XTE J1550-564, Miller et al. (2001); 240, 165 Hz in 
H1743-322, Homan et al. (2005)) were used to estimate the parameters of the BHs such as 
the masses and the spins. In fact, the frequencies are stable, i.e. there are small varitions in 
them, rather than invariable (Miller et al. 2001; Remillard et al. 1999; Morgan et al. 1997; 
Strohmayer 2001b). 

This paper is organized as follows. In Section 2 we suggest the basic model and get two 
stable modes of the Alfven wave by solving the general relativistic magnetohydrodynamic 
(GRMHD) equations of the perturbed plasma in BH accretion disks. In Section 3 we compare 
the results with the observations and discuss their possible implications. 

2. The GRMHD model of the HFQPOs in BH-LMXBs 

We consider that the HFQPOs of BH-LMXBs result from GRMHD waves caused by 
perturbations in the disk. According to §2.1 we can find that only one type of GRMHD wave, 
the steady Alfven wave, can spread along the magnetic field lines to the energy release region. 
The simulation by Koide (2002) also shows that a torsional Alfven wave can be generated 
by the rotational dragging of space. We assume that the puny disturbance doesn't change 
the metric of the space-time. 
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2.1. The modes of the GRMHD wave 

The oscillation modes of plasma in accretion disks with or without magnetic field, such as 
g- modes, p-modes, and c-modes, have ever been investigated (Wagoner et al. 2001; Fu & Lai 
2009; Lai 2009). Here we discuss the GRMHD wave modes in ideal adiabatic magnetofluid. 
In the fiducial observer (FIDO) frame, which is a locally inertial frame, the line element can 

3 

be written as (ds) 2 = — (cdt) 2 + ^ (da;*) 2 , where c is the speed of light in vacuum, the Roman 

i=i 

indices (i) run from 1 to 3, and (ct, x 1 , x 2 , x 3 ) are the coordinates of the FIDO frame. The 
metric of the Kerr space-time in the Boyer-Lindquist coordinates (ct' ,r,9,ip) is introduced 
and the line element for the observer at infinity (i.e. in the laboratory frame) can be written 

as, 

3 

(ds) 2 = -hl(cdt'f + [hl^'f ~ 2h 2 0Oidt'dx% (1) 



i=i 



where 



hi = 1 - ^,hl = |, hi = £,/> 2 = ^sin 2 9, u 1 = u 2 = 0, c 3 = (2) 

Here (ct', x 1 ', x 2 ', x 3 ') are the coordinates of the laboratory frame, ho, hi and Ui are the 
metrics of the Kerr space-time, r g = GM/c 2 , a = Jc/GM 2 (M and J are the mass and the 
angular momentum of the BH respectively), G is gravitational constant, £ = r 2 + a 2 r 2 cos 2 9, 
A = r 2 — 2rr g + a 2 r 2 , and A = (r 2 + a 2 r 2 ) 2 — a 2 r 2 sin 2 9 A. Since our discussion is limited 
to the accretion disk we take 9 = n/2 and get S = r 2 and A = (r 2 + a 2 r 2 ) 2 — a 2 r 2 A, where 
r is the distance of the plasma from the BH. Besides that, the lapse function and the shift 
velocity can be expressed as, 



a 



^ Ac^ / r 4 - 2r% + a 2 r 2 r 2 

,l "T V r ) \ 4 _|_ , O n 2 rr 3 ' 



c V r 4 + a 2 r 2 r 2 + 2a 2 rr^ ' 



and 



or /3 = (A, /3 2 , ft) = (0, 0, g == ), (4) 

ca r y/r z — 2rr 9 + a^r^ 

where /3 is a vector parallel to the toroidal velocity of the plasma. 

We begin with the form of 3+1 split of the GRMHD equations as in Koide (2003). 
Previous investigations (Ruzmaikin et al. 1979; Tout et al. 1992; Ruediger et al. 1995; 
Hawley et al. 2000; Moss et al. 2004; Hirose et al. 2004) have shown that in the accretion 
disk around a BH the toroidal component of the magnetic field may be much stronger than 
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the poloidal component, i.e., B r <C B v and B$ <C B v . In our analysis we also assume 
v r <C and vg <C v v . The equations in the FIDO frame are then written as follows: 

^ = -V> 7 p(v + c/3)], (5) 
<9P ~ 

— = -V • [a(T + c/3P)] - (e + 7 pc 2 )Va + af curv - P • a, (6) 
at 

de ~ 

= -V • [«(c 2 P - 7 pc 2 v + ec/3)] - (Va) • c 2 P - T : a, (7) 

— = -Vx [a(E-c/3xB)], (8) 

V • B = 0, (9) 

E + vxB = 0, (10) 

pp~ T = constant, (11) 

where v is the velocity of the plasma, p the plasma density, p the barometric pressure, T 
the adiabatic index, 7 the Lorentz factor, E = E'/\/pO; B = B'/\/A*o (here E' is the electric 
field, B' the magnetic field, and pio the magnetic permeability in the vacuum), respectively. 
The bold characters denote vectors, and the superscript ~ corresponds to tensors. The 
energy-momentum tensor is, 

T = (p+^ + ^)I+ | 7 2 VV - BB - 1 EE , (12) 

where ip = pc 2 + ^Py is the relativistic enthalpy density. The equivalent momentum density 
and energy density are, 

P = | 7 2 V + ^E x B. (13) 

and 

B 2 E 2 

e = W -p-ipc 2 + — + — , (14) 

respectively. Equations (5), (6), (7), (8), and (9) correspond to the continuity equation, the 
momentum conservation equation, the energy conservation equation, the law of electromag- 
netic induction, and the equation of no divergence, respectively. Equations (10) and (11) 
are the equations of state for the infinite electrical conductivity and the adiabatic condition. 
The other physical quantities are, 

f curv = .( G ij Tij ~ GjiT jj ), a l3 ee (hi/h^doui/dx^, Gn = -(1/hihj) x (dh t /dx>). (15) 
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We can simplify the above two tensors in the accretion disk in the Kerr space-time as, 








2a r 

2acr 2 (3r 2 +a 2 r 2 )]l r 2 +a 2 r 2 + 










^ + a 2 r 2 -2rr g (r 3 +ra 2 r 2 +2a 2 rl) 2 




°). ( 16 ) 



and 



G, 



r 2 r 2 — 2rr g +a 2 r 
_ J y/r 2 -2rr g +a 2 r 2 



(r 3 —a 2 r'g) ^Jr 2 —2rr g +a 2 r 2 
(r 3 +ra 2 r 2 +2a 2 r' g v )r 2 










(17) 



where we define the vector N = [(h 3 /hi)(du 3 /dx 1 ), 0, 0]. The magnetized accretion torus 
with a toroidal magnetic field around a Kerr black hole can exist stably (Komissarov 2006), 
and the GRMHD equations in the steady state can be expressed as, 



d(7Po) 
dt 



-V • [«7p (vo + c/3)], 



<9P 

dt 

dt 



-V • [a(T + c/3P )] - (e + lp c 2 )Wa + af curv , - P • <r, 

: -V • [a(c 2 P - 7Poc 2 v + e c/3)] - (Va) • c 2 P -T :a, 
dB 



dt 



-V x [a(E -c(3 x B )], 

V • B = 0, 
E + v x B = 0, 

PoPo r = constant, 



(18) 

(19) 

(20) 

(21) 
(22) 
(23) 
(24) 



where the subscript denotes the variables in steady state. Next we consider the GRMHD 
equations after the plasma is perturbed slightly, 



dt 



= -V- [a 7 p(v + c/3)], 



dP 

~dt 
d£ 
dt 



-V • [a(T + c/3P)] - (e + 7 pc 2 )Va + af curv - P • a, 
-V • [a(c 2 P - 7pc 2 v + ec(3)} - (Va) • c 2 P - f : a, 



(25) 

(26) 
(27) 
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— = -Vx [a(E-c/3xB)], (28) 

V • B = 0, (29) 
E + vxB = 0, (30) 

PP~ T = PoPo r ? ( 31 ) 

where v = v + v s , p = p + p s , p = p + p s , e = e + e s , E = E + E s , B = B + B s , 

P = Po + Ps, /curv = fcurvfl + jcurv.s, and T = T + T s , with the subscript s denoting the 
perturbed quantities (v s <C v , E s <C E , B s <C B Q , p s <C p Q , and p s <C Po) and with the 
superscript " denoting the variables after the disturbance. By combining Eqs. (18)-(31) we 
get the equations about the perturbed quantities in the first-order approximation, 

d{lPs) 



dt 

dP 



- V • [a 7 (p v s + p s v + p s c(3)] , (32) 



s V • [a(T s + C/3PJ] - (e s + 1Ps c 2 )Va + a/ curv , s - P s • a, (33) 



dt 



de 

-± = _ V ■ [a{c 2 P s - 7 c 2 p v s - 7c 2 p s v + e s c(3)} - (Va) • c 2 P s -T s :a, (34) 

-^ = -V x[a(E s -cf3xB s )], (35) 

V-B s = 0, (36) 
E s + v xB s + v s xB„ = 0, (37) 

p s = p s . (38) 

Po 

With Eqs. (12)-(14) we can get the perturbed energy-momentum tensor (T s ), equivalent 
momentum density (P s ) and equivalent energy density (e s ) in the same way, by defining 
ifo = Poc 2 + fzi and ip s = p s c 2 + 

T s = (p s + B • B s + 1e • E S )I + ^7 2 V V + ^7 2 (V V S + V S V ) 

-(B B s + B s B )-^-(E E s + E s Eo), (39) 
& 

7 2 1 1 

P s = -^(^ov s + ^ s v ) + -E s x B + ^E x B s , (40) 
e = ^ sl 2 -p s - 1Ps c 2 + B • B s + (41) 
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After carrying out Fourier transformation (e lk ' r luJt ) for Eqs. (32)- (36) and substituting 
Eqs. (37) and (38) into them, we get the following equations when v ||B ||/3 is considered, 

[u - a(k ■ v ) - ac(k ■ (3)}p s = ap (k ■ v s ), (42) 

ujP s - ac(k ■ (3)P S = ak ■ T s + a(e s + 7p s c 2 )k + iaf CUTVtS - iP s ■ a, (43) 

ue s - ac(k ■ (3)e 8 = 2ac 2 (k • P s ) - a>yc 2 p (k ■ v s ) - a^c 2 ik • v )p s - iT s : a (44) 

[uj - a(k ■ v ) - acik ■ f3)]B s = a(k ■ v s )B - a(k • B )v s , (45) 

k • B s = 0, (46) 

where k is the wave vector and oo is the the oscillation frequency. When Eqs. (37) and 
(38) are substituted into Eqs. (39)-(41) and v ||B ||/3 is considered, Eqs. (39)-(41) can be 
converted to be, 

T s = {^p s + B • B S )I + ^7 2 V V + ^7 2 (V V S + V S V ) - (B B S + B S B ), (47) 
Po c c 

P s = ^(^ov s + ^ s v ) - l(v • B )B S + 1(B S • B )v - ^(v s • B )B + (48) 

c c c c c 

e = Vs7 2 - — Ps - ip s c 2 + B • B s . (49) 
Po 

If lu = ak ■ v + ack ■ {3, we get a unphysical solution (lu = 0) from Eqs. (42)- (49), because 
the conditions k • v s = and k • v = are derived from Eqs. (42)-(44). Now we consider 
w^ak-v + ack ■ (3 and discuss all the three types of MHD waves: the Alfven wave, the 
ion-acoustic wave, and the magnetosonic wave. 



2.1.1. The Alfven wave 

The transportation direction of Alfven waves is along the magnetic field line, so kff v ||B ||/3, 
where || denotes parallel and || the same direction. Since the Alfven wave is a transverse 
wave, i.e. k _L B s , k _L v s or the ip component of the perturbed velocity, v Sj(p = 0. 

From Eqs. (15)-(17), (42), (47), and (48) and the above results we obtain 

fcurv, S = (0,0,0), 
7 2 

P s -cr = -^WoV a , v + VwN = 0, 

and 

^ ~ r 7 2 , otkBl 
T : a = a 3 i[—^ v + 



c 2 uj — aik • v ) — ac{k • (3) 
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where v Sjr denotes the r component of the perturbed velocity. Substituting the above equa- 
tions and Eqs. (42), (45)- (49) into Eqs. (43) and (44) we can get v Sjr = and the dispersion 
equation, 

[w-ac(k-/3)]% 2 ^o + £o)-( 2 a7^o^ = 0. (50) 

The frequencies of the Alfven waves are solved as, 



u = ka[^c + WT - Wo ]. (51) 

The group velocities of the Alfven waves are in the same form of the phase velocities of the 
Alfven waves and are, 



uj _du _ l 2 j>oVo ± B^B 2 c 2 + (c 2 - gggo n ( - , 

k~dk~ a[hC + Bl + 7^ J - (52) 

These velocities can be simplified in the special relativity, i.e. r — >■ oo, f3 3 — >■ & a — >■ 1, as 



^° ± 71 + (f) 2c _ Uo ± ^V?+^ C 



where r y 2 r] 2 = B^/ip®. Equation (53) is in the same form as the expression of De Villiers & 
Hawley (2003). 



2.1.2. The ion-acoustic wave and the magnetosonic wave 

The ion-acoustic wave is a longitudinal wave without electromagnetic polarization so 
k||v s and B s = 0. According to Eq. (45) we have Bo||V s owing to B s = 0. This leads to the 
conclusion that k||v s ||B ||v ||/3, and that no ion-acoustic wave solution fits Eqs. (42)-(49). 

The magnetosonic wave is another type of longitudinal wave with transverse electromag- 
netic polarization, i.e., k || v s , k _L B s . The wave vector k isn't parallel with B under the 
condition Bo||vo||/3, otherwise the wave is degenerated into the ion-acoustic wave and does 
not exist. Now we discuss the magnetosonic wave under two different conditions respectively. 

Firstly when k _L v we get k || v s , k _L B s , B \\v \\(3 _L k, and Eqs. (42)-(49) can be 
simplified to be, 

(«7 2 ^o + cxTp Q - cvy 2 Tp + <xB 2 )(k ■ v s ) = iT s : a, (54) 
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r 2 

w 2 (7 2 V ; o+5o)v s +7 2 arp (k-v s )wv -« 2 c 2 /c 2 (25o+7 2 poC 2 +7 2 ^-^po)v s = ic 2 uj(ai CUIV:S -P s -a), 

(55) 

~ _ 7 2 h 3 du 3 

T ° :a = ^h^ mVs ^ (56) 

P s • of = ^a( Po c 2 + T ArPo)-(k • v s )N, (57) 
c z 1 — 1 u; 

where 

7 2 

fcurv.s = (fs,v, 0, G31— ^0^0^,r), (58) 

with 

/ s ,r = -G 2 i«(k • v s ) ^ - G 31 [a ± + -5- a(po<r + - — -p )— ](k • v s ). 

10 U! C 1 — 1 U! 

From Eqs. (54)- (58) we can obtain an unstable solution (an increasing wave or an attenuation 
wave) if 

T s : a c 2 f s ^ k • v s 



«7 2 ^o + «r> - «7 2 rpo + oiBl r7 2 p fo 
or no solution when k _L v . The possible unstable solution is 



A + A 2 + 4( 7 2 ^o + Bl)(2Bl + 7 2 p c 2 + f^- Po )a 2 c 2 k 

2 7 2 ^ + 2B 2 

where 



u = 1 ■ 002 . ( 59 ) 



T^poVoVs k • v s 

and (P s • cr) r represents the r component of the vector P s ■ a. 

Secondly when k is neither parallel with nor vertical to v we get the solution from 
Eq. (42)-(49) as follows, 

za(k-v )(T s :3 : ) + 2a 2 7 2 rpo(k-v s )(k-vo) 
uj = ac(k ■ p) - , 

(a 7 2 ^o + al> - «7 2 rp + aB 2 ) (k • v s ) - a(k • B ) (B • v s ) -iT s :a 

(60) 

which is also unstable. 

We summarize the results in §2.1.1 and §2.1.2: (1) there are two stable Alfven wave 
modes, (2) the ion-acoustic wave doesn't exist, and (3) a few unstable modes of magnetoa- 
coustic wave may emerge in the GRMHD in the BH accretion disks. 
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2.2. The relation between the magnetic energy density and the relativistic 

enthalpy density 

It is widely believed that the toroidal magnetic field in accretion disks is generated by 
dynamo mechanism (Ruzmaikin et al. 1979; Tout et al. 1992; Ruediger et al. 1995; Hawley 
et al. 2000; Moss et al. 2004), and that accretion is driven by the magnetic stress (e.g. 
Matsumoto & Tajima 1995; Brandenburg et al. 1995; Stone et al. 1996). Accordingly the 
angular momentum conservation gives (Torkelsson 1998) 

B' B' 

Mr Vip = (61) 
Po 

where H is the thickness of the accretion disk, B^ and B' r are the </?- and r-components 
of the magnetic field B', M(= 27iHrv r p) is the accretion rate, v v and v r are the ip- and 
r-components of the velocity of the accreting plasma, respectively. Suppose B'^ = 7d yn -5r 
and v r = lv v , the above equation can be simplified to be, 

% 2 = ^P/^oTdyn- ( 62 ) 

The velocity v v can be approximatively expressed as the velocity of the circular orbit relative 
to Bardeen observers, which can be got from the equations (8.354-8.359) of Camenzind 
(2007). The velocity of the plasm in the prograde orbit is, 

v v = [(i 2 + a 2 + -i a l) i3 ^ 2i 2^ i+a ti^ i+ati 2 +a 3^ 2afi - f a *] v / i 2_ C 2i+a 2' ( 63 ) 
and in the retrograde orbit, 

v = [( ? 2 + al + la*) . 3 T . ^f~fr al ~T s „ ■ - -«*] / ° , (64) 

where a* = |a|, % = r/r g . The two velocities return to the Keplerian velocity c/y/i when 
i — ?■ oo. 



Since frypo <^ Poc 2 in the accretion disks, we can express = 7 2 ?7 2 (poc 2 + ^zj) as 

b' 2 

= _B 2 ~ 7 2 ^ 2 poC 2 and ^ 2 rf = 7d yn ^ 2 /c 2 = l'v 2 /c 2 from Eq. (62), where I' = 7 dyn / = 

B 1 

-gf We then get the relation between the magnetic energy density and the relativistic 



enthalpy density, ^ ~ ^■f^^rV'o, and express the frequencies of the Alfven waves, i.e., 
Eq. (51), as, 

7 2 ± Jr + r 2 ! 

u = ka[P 3 c H -2 v v \. (65) 

l 2 + l'^ 
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When we consider the characteristic wavelength A ~ r , i.e. the wave number k ~ 27r/r, 
Eq. (65) is reduced to be 

u ~ — c£ = — c£ (66) 



r «r 9 



where £ = a p 3 H * - ,, 2 . \ * 



2.3. Estimate of the parameter /' 

Begelman and Pringle (2007) have investigated the structure of accretion disks with 
strong toroidal magnetic fields, and found that the thickness of the disks is higher than 
standard thin disks, but in line with observations (Robinson et al. 1999; Shatter & Misselt 
2006), 

| = 0.48a-/i 7r ^Mr (^)" 15/68 , (67) 

where a* is the viscosity prescription, M 18 = M/10 18 gs _1 , and r w = r/10 10 cm. The 
corresponding radial velocity is, 

3 H 

Vr ^ -a*(—) 2 v k , (68) 
2 r 

where is the Keplerian velocity. Here we adopt that v v ~ to estimate the value of I'. 
From Eqs. (67) and (68), the parameter / is, 

1=** 0.3456«f/ 17 rr M^ 17 (^)- 15/34 . (69) 

Vishniac et al.(1990) have discussed the dynamo action by internal waves in accretion disks 
and suggested B T /B V ~ a*, so we have 

I' = / 7dyn - 0.3456«; 2 / 17 rr M^ 17 (^)" 15 / 34 . (70) 

The relatively small values of the power indices in Eq. (70) indicate that I' is not sensitively 
dependent on a*, r, M and M. 

HFQPOs in BH-LMXBs were generally observed in the steep power-law (SPL) state, i.e. 
very high state (VHS), and it is very likely that the QPOs are associated with the region that 
produces hard X-ray emission. The accretion model for the VHS is a long-debated subject. 
It might contain an inner ADAF surrounded by an outer thin disk (e.g. Yuan 2001). By 
analyzing the observational data of GRO J1655-40 and XTE J1550-564, McClintock & 
Remillard (2006) showed that the disks with blackbody radiation appear to truncate at a 
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radius (~ 100r g ) in the low/hard state, and that the truncated radius decreases when the 
power-law component becomes stronger and steeper. Theoretical investigations (Abramowicz 
et al. 1995; Honma 1996; Esin et al. 1997; Liu et al. 1999; Manmoto et al. 2000; RozaNska 
et al. 2000; Narayan et al. 2008) also suggest that the transition radius ranges from ~ 100r g 
to 10000r g . Accordingly the transition radius in VHS is likely to be smaller than in the 
low/hard state, and for illustration here we adopt its value to be 50r g . The X-ray luminosity 
in the VHS is often more than 0.2L Edd where L Edd is the Eddington luminosity (McClintock 
& Remillard 2006), so we adopt a typical accretion rate 10 18 gs _1 . We also take a* ~ 0.1 
(King et al. 2007). From the above values we can estimate the parameter /' as, 

/' „ ,048(^-)- 2 / 17 (^-) 9 / 34 ( 1R M )3/17(JL)-15/34 (n) 

V 0.1 ; v 50r g V 10 18 gs" i; V 7M ; 
Combing Eqs. (3), (4), (66), and (71) we get the two frequencies of the HFQPOs as, 

w = 1.2756 x lO 6 ^^-)- 1 (72) 

1 Mq 



where 



2*3 + M r 2a 7^0.04765 + 0.0022705^/^ 



i 4 + a 2 i 2 + 2a 2 i L iV« 2 - 2i + a 2 7 2 + 0.04765^/c 2 c J ' 

Our calculations show that, for BH-LMXBs with measured HFQPOs the ratio of the two 
frequencies is generally around 1.5 (see Table 1). 



3. Results and discussion 

Equations (63), (64), and (72) indicate the existence of Alfven waves with two fre- 
quencies in the accretion disk. The ratio of the upper and lower frequencies is close to 
3 : 2, suggesting that these waves may account for the HFQPO pairs. Given the mass and 
spin of a BH, one can determine the radius where the QPOs are produced. In Table 1 we 
present the inferred radius by comparing Eq. (72) with either the upper or the lower centroid 
QPO frequency of several BH-LMXBs. We adopt the averaged masses for GRO J1655— 40 
(6.0 - 6.6M , McClintock & Remillard 2006), GRS 1915+105 (10 - 18M , Greiner et al. 
2001), and XTE J1550-564 (8.4 - 1O.8M , McClintock & Remillard 2006) and averaged 
dimensionless spins of GRO J1655-40 (0.65 - 0.75, Shafee et al. 2006) and GRS 1915+105 
(0.98 - 1.0, McClintock k Shafee et al. 2006). Because the spin of XTE J1550-564 hasn't 
been measured, we take it to be 1, 0.5, and 0. It is interesting to see that in most cases the 
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radii are ~ 70r g , consistent with the transition radii between an ADAF and a thin disk dis- 
cussed above. The frequencies of the HFQPOs mainly depends on the transition radius (r tr ) 
and the ratio of I' = (v I /v (p )/(B' r /B' ). According to the discussion in §2.3 the parameter /' 
changes little with the a*, r, M, and M, this may explain why the HFQPO frequencies are 
relatively stable during the VHS. 

Considering the similarities in BH and NS accretion disks, it is interesting to ask why 
the 1.5 frequency ratio is not evident in the HFQPOs in NS-LMXBs. The reasons may lie in 
the differences in the configuration of the magnetic fields and the structure of the accretion 
disks in BH- and NS-LMXBs. It is well known that NSs generally hold dipolar magnetic 
fields, and the toroidal field component, induced by difference in the angular velocity between 
the disk and the NS magnetosphere, can never be stronger than the poloidal component, 
otherwise the field configuration becomes unstable leading to inflation of the field lines (e.g. 
Aly 1985). On the other hand, both theories and simulations (Ruzmaikin et al., 1979; Tout 
et al., 1992; Ruediger et al., 1995; Hawley et al., 2000; Moss et al, 2004; Hirose et al., 2004) 
show that in BH accretion disks the toroidal magnetic field can be much stronger than the 
poloidal one. Additionally, recent observations suggest that the inner disk around an NS 
may not be an ADAF even in low state, as the thermal emission from the surface of the NS 
would tend to cool an ADAF, making such a flow more difficult than in BH systems (Cackett 
et al., 2009). 

Finally, a correlation between the lower HFQPO frequencies and the BH masses has 
been suggested, i.e., v\ oc M~ l (McClintock & Remillard 2006). This correlation can be 
naturally reproduced in our model. In Fig. 1 we plot the predicted relation between v\ and 
m when r tr changes from 66r g to 76r g , which fit reasonably with the measured data. 
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by the Natural Science Foundation of China (under grant number 10873008) and the National 
Basic Research Program of China (973 Program 2009CB824800). 
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Table 1: This table presents the masses (m in solar units), dimensionless spins a for several 
BHs, and the inferred radii (r tr u and r tr l ), the predicted frequency ratios (zAi/Vi) in the 
corresponding sources based on the measured upper and lower QPO frequencies (z/ u and v\ 
in Hz) respectively. 



sources 


m 


a 




r t r,u 






r tr ,i 




GRO J1655-40 


6.3 


0.7 


450 


66.7528 


1.5504 


300 


65.3002 


1.5501 


GRS 1915+105 


14 


0.99 


168 


75.6114 


1.5514 


113 


73.4984 


1.5510 


GRS 1915+105 


14 


0.99 


67 


139.7920 


1.5570 


41 


144.368 


1.5572 


XTE J1550-564 


9.6 


1 


276 


69.8141 


1.5503 


184 


68.2975 


1.5500 


XTE J1550-564 


9.6 


0.5 


276 


69.8639 


1.5515 


184 


68.3125 


1.5512 


XTE J1550-564 


9.6 





276 


69.9159 


1.5526 


184 


68.3302 


1.5524 
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Fig. 1. — The relation between the lower HFQPO frequencies and BH masses for GRO 
J1655-40, XTE J1550-564, and GRS 1915+105. 



